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Comprehensive x-ray scattering studies have characterized the smectic ordering of octyl- 
cyanobiphenyl (8CB) confined in the hydrogen-bonded silica gels formed by aerosil dispersions. 
For all densities of aerosil and all measurement temperatures, the correlations remain short range, 
demonstrating that the disorder imposed by the gels destroys the nematic (N) to smectic-A (SmA) 
transition. The smectic correlation function contains two distinct contributions. The first has a form 
identical to that describing the critical thermal fluctuations in pure 8CB near the N-SmA transition, 
and this term displays a temperature dependence at high temperatures similar to that of the pure 
liquid crystal. The second term, which is negligible at high temperatures but dominates at low 
temperatures, has a shape given by the thermal term squared and describes the static fluctuations 
due to random fields induced by confinement in the gel. The correlation lengths appearing in the 
thermal and disorder terms are the same and show strong variation with gel density at low tempera- 
tures. The temperature dependence of the amplitude of the static fluctuations further suggests that 
nematic susceptibility become suppressed with increasing quenched disorder. The results overall are 
well described by a mapping of the liquid crystal-aerosil system into a three dimensional XY model 
in a random field with disorder strength varying linearly with the aerosil density. 



I. INTRODUCTION 

Liquid crystals have long served as important model 
systems in statistical mechanics. For example, experi- 
ments on phase transitions in liquid crystals have pro- 
vided many of the most detailed tests of the modern 
theories of critical phenomena. Recent studies in liq- 
uid crystals have investigated the effects that quenched 
disorder produces in phase behavior and mesophase or- 
dering, and a fruitful strategy in experiments for intro- 
ducing quenched disorder has been through confinement 
in random porous media. The fragile nature of the me- 
somorphic phases and the importance of surface interac- 
tions make such confinement particularly well suited for 
liquid crystalline systems. For example, a series of exper- 
imentsETEl has characterized the effects on the nematic to 
isotropic and nematic to smectic-A transitions of liquid 
crystals confined in aerogels - highly-porous, chemically 
bonded silica gels. The primary conclusion of this work 
has been that the transition behavior becomes severely 
smeared and that long-range nematic and smectic order- 
ing is suppressed by the aerogel. In particular, these 
studies have verified theoretical expectations about the 
fragility of the smectic phase to disorden3. 

In order to access a weaker regime of disorder than is 
possible with aerogel or other rigid porous media, several 
recent studies have focused on liquid crystals confined 



within aerosil gels - weak, thixotropic gels comprised of 
nanometer scale silica particles. The disorder imposed 
by the aerosil gels can be made less severe than in the 
case of aerogels both because lower volume fractions are 
possible with aerosils and because the compliance of the 
aerosil gels leads to partial annealing of the disorder. 
This work on liquid crysjtals in aerosil dispersions, includ- 
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dielectric susceptibility^, and 
has revealed important differ- 
ences between such systems and liquid crystals in aero- 
gels, which can be ascribed to the weaker nature of the 
aerosil disorder. In particular, an x-ray scattering inves- 
tigation into the nematic to smectic-A transition in a pro- 
totypical thermotropic liquid crystal, octylcyanobiphcnyl 
(8CB), in the presence of dispersed aerosil gels has shown 
that the aerosil gel can be quantitatively understood as 
introducing weak-to- intermediate strength random fields 
that, compete with smectic ordering in the liquid crys- 
tafi 

Random field sytems have been fruitful models for ex- 
ploring effects of quenched disorder, and experimentally 
random field Ising magnets have been at the center of 
this workE3. The nematic to smectic transition breaks a 
3D-XY symmetry, thus 8CB with dispersed aerosils pro- 
vides an opportunity to study experimentally ajandom 
field system that breaks a continuous symmetryEjlia. In 
this paper we provide a comprehensive picture of the 
smectic correlations in 8CB confined in aerosil gels, as 
determined by the x-ray scattering. The analysis we ap- 
ply adopts a random-field picture, and consistent with 
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theoretical predictions for a 3D-XY system with random 
fields, we find that the smectic phase is destroyed by 
the disorder and is replaced by the growth of short-range 
smectic correlations. A noteworthy feature ofrthese x-ray 
results and corresponding calorimetric studiesQ is the lack 
of any measurable hysteresis or time dependent effects 
that would indicate out-of-equilibrium behavior. In the 
random field Ising magnet slow dynamics and metasta- 
bilityoEj severely complicate efforts to understand the 
underlying equilibrium behavior. The absence of such 
problems for the smectic in aerosil, presumably a conse- 
quence of the system's continuous symmetry, thus pro- 
vides a unique perspective on equilibrium behavior in a 
random field system. The companion paper that follows 
this papertll, hereafter called Paper II, compares the x- 
ray and calorimetry results to reveal scaling behavior, 
including finite-size scaling effects, in the smectic order- 
ing due to the presence of the aerosil gel. 

Recent theoretical work has emphasized the strong 
effect of confinement in random media on liquid crys- 
tal phases^ 'BErO. A detailed study by Radzihovsky 
and TonertH, which models the confinement as introduc- 
ing random fields, has predicted the destruction of the 
smectic-A phase by arbitrarily weak quenched disorder, 
consistent with experiment, and has introduced the pos- 
sibility that a topologically ordered "smectic Bragg glass" 
phase may appear at low temperatures. The smectic 
Bragg glass is a possible manifestation of the anoma- 
lous elasticity that the theory predicts for the smectic in 
the presence of disorder. Distinguishing features of this 
anomalous elasticity and smectic Bragg glass phase are 
the smectic correlations in the system. Thus, x-ray scat- 
tering, which can directly probe the smectic correlation 
function, is an ideal probe for testing some of these ideas, 
and recent analysis of x-ray studies on 8CB in aerogels 
has been reported to show agreement with many of these 
predictionscl. The weaker nature of disorder from aerosil 
gels should place 8CB + aerosil samples more firmly in 
the regime addressed by the theory. However, as shown 
in the present work, the agreement between the observed 
smectic correlations in 8CB in aerosil gels and more de- 
tailed predictions of the theory is limited. 

Section II of this paper describes the sample prepa- 
ration and the details of the x-ray scattering measure- 
ments on 8CB in aerosil gels, and Sec. Ill outlines the 
results of the x-ray lineshape analysis, which yields in- 
sight into both thermal and static smectic fluctuations 
and the temperature dependence of correlation lengths 
and smectic susceptibilities. Section IV provides a dis- 
cussion of the pseudo-critical structural behavior of 8CB 
in weak aerosil gels. The latter aspects are extended by 
the scaling behavior and comparisons with calorimetry 
covered in Paper II. 



II. EXPERIMENTAL PROCEDURES 
A. Sample preparation and characteristics 

The 8CB used in this study came from a single syn- 
thetic batch obtained from Aldrich Corp. The liquid 
crystal had a quoted purity of 99% and was used with- 
out further purification. In the absence of disorder, pure 
8CB undergoes an isotropic (I) to nematic (N) transition 
at Tjlfj = 313.98 K and a nematic to smectic-A (SmA) 
transition at T% A = 306.97 Ml Below 290 K, 8CB forms 
a three-dimensional ordered crystal phase. The aerosil, 
obtained from Degussa Corp., consists of 70-A diameter 
SiC>2 spheres, and the type 300 aerosil used in this study 
is strongly hydrophilic. 

The aerosil gels were formed directly in the liquid crys- 
tal following established procedures^. Appropriate quan- 
tities of degassed 8CB and dried aerosil powder were 
mixed with high purity acetone, and the suspension was 
sonicated for several hours to achieve a uniform disper- 
sion. The suspension was then gently heated to 315 K to 
evaporate the acetone slowly. After no signs of acetone 
remained, the samples were placed under vacuum at 10 -2 
Torr for 12 hours at 320 K to remove any trace amounts 
of solvent or absorbed water vapor. For aerosil densities 
above a gelation threshold of approximately 1% silica by 
volume, the resulting material was a highly uniform soft 
solid that maintained its shape when heated above T^j. 
Small-angle x-ray scattering has revealed that aerosil gels 
formed under trus-.procedure bear a strong resemblance 
to rigid aerogelaJ'Q. An extended discussion of the gel 
structure formed by aerosils is given in Paper II. The 
samples in this study ranged in aerosil densities from 
ps = 0.025 g sil/cm 3 8CB, which is just above the gela- 
tion threshold, to ps = 0.341 g sil/cm 3 8CB. Several 
efforts to prepare samples below the threshold density 
produced macroscopically inhomogeneous materials that 
were unsuitable for study. 

B. X-ray scattering 

The x-ray scattering studies were conducted on the 
X20A and X20C beamlines of the National Synchrotron 
Light Source at Brookhaven National Laboratory using 8 
keV x-rays. The samples were placed in aluminum hold- 
ers with epoxy-sealed Kapton windows designed to main- 
tain a sample thickness of approximately 1 mm, closely 
matching the attenuation length for the 8 keV radiation 
through the material. The holders were mounted in a 
brass block which in turn was positioned in a beryllium 
can containing dry nitrogen gas. A thermoelectric cooler 
mounted on top of the beryllium can allowed us to cool 
the sample below room temperature when desired, and a 
home-built P-I temperature controller maintained a set 
temperature for the brass block. Preliminary x-ray scat- 
tering measurements showed that the crystallization of 
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FIG. 1. Scattering intensity for 8CB in an aerosil gel with 
ps = 0.161 g/cm 3 at three temperatures spanning the N-SmA 
transition of the pure system, T% A — 306.97 K. The curves 
for 300.7 K and 306.7 K have been shifted by factors of 100 
and 10, respectively, for clarity. The strongly temperature 
dependent peak near q = 0.2 A" 1 corresponds to scattering 
from smectic fluctuations. 

8CB was suppressed approximately 7 K by confinement 
in aerosil gels. Because such crystallization irreversibly 
damages the aerosil gels, we limited our studies to tem- 
peratures safely above this point. The measurement tem- 
peratures ranged from 287 K, approximately 20 K below 
the N-SmA transition in pure 8CB, to 318 K, approxi- 
mately 4 K above the N-I transition. The temperature 
stability during a measurement was better than ± 0.001 
K. 

The scattering intensity was measured in transmission 
in a vertical scattering geometry. The scattered beam 
was reflected from a single-crystal germanium analyzer 
to achieve high wavevector resolution and was measured 
with a point scintillation detector. Because the scatter- 
ing intensity from the liquid crystal-aerosil gel compos- 
ites had the azimuthal symmetry of a powder pattern, 
two sets of slits between the sample and the analyzer 
were used to define tightly the horizontal acceptance and 
thus provide an undistorted measurement of the scatter- 
ing intensity as a function of wavevector transfer, I(q)cB. 

Preliminary measurements of the scattering from 8CB 
+ aerosil samples revealed damage to the material from 
prolonged exposure to the x-ray radiation. This dam- 
age was reflected in irreproducibility in the observed I{q) 
and a progressive shift in the onset of the short-ranged 
smectic order described in Section HID. To combat this 
problem a protocol was developed which made these ef- 
fects immeasurably small. Immediately after the sample 
was loaded into the scattering cell and before any x-ray 
measurements were made, the material was heated to 
approximately 318 K, several degrees above T^j. The 
material was then cooled through a series of tempera- 
tures at which measurements of the scattering intensity, 
I(q), were made. A process of trial-and-error during pre- 
liminary studies revealed the appropriate exposure times 
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FIG. 2. Comparison of the smectic Bragg peak measured 
for 8CB in aerosil gels with that measured for pure 8CB. All 
peak heights have been normalized to 1. The solid line is for 
pure 8CB below the N-SmA transition temperature. The solid 
circles correspond to ps = 0.025 g/cm 3 and the open circles 
to ps = 0.220 g/cm 3 . The width of the peak for pure 8CB 
is essentially identical to that set by the experimental resolu- 
tion, as determined from the direct x-ray beam profile. The 
wider peaks measured for 8CB in aerosil gels indicate that the 
smectic correlations remain short-ranged to low temperatures 
(approximately T% A — 19 K shown here) due to the disorder 
imposed by the gels. 

for the final measurements to optimize statistics while 
avoiding measurable damage. After the series of mea- 
surements ending at the lowest temperature was com- 
plete, the material was reheated and measurements were 
repeated at several temperatures to assure that no mea- 
surable changes in the scattering had occurred. Finally, 
the material was heated back to 318 K, and a measure- 
ment with a long counting time was conducted to de- 
termine the background scattering from the aerosil gel 
structure. 



III. RESULTS 
A. Lineshape analysis 

Figure 1 displays representative results for I(q) from 
a sample with ps = 0.161 g/cm 3 at three tempera- 
tures spanning the N-SmA transition temperature in 
pure 8CB, T% A = 306.97 K. The scattering lineshape 
contains two salient features: a broad, sloping back- 
ground that is approximately temperature independent 
and a strongly temperature-dependent peak that devel- 
ops on cooling near 0.2 A" 1 . The background consists 
primarily of scattering from the aerosil gel structure. The 
temperature-dependent peak results from smectic fluctu- 
ations in the 8CB. The shape of this peak is a direct 
measure of the smectic correlation function for the 8CB 
confined in the aerosil gel. 

For the full range of aerosil densities and to the lowest 
measurement temperatures, the observed smectic peak 
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is significantly broader than the resolution limit set by 
the x-ray optics, as illustrated in Figure 2. Measure- 
ments of pure 8CB below T^ A , shown by the solid line 
in the figure, revealed a lineshape whose width is virtu- 
ally indistinguishable from that of the resolution, which 
we determined from the shape of the profile of the direct 
x-ray beam. (The subtle increase in width in pure 8CB 
due to Landau-Pieirls instability is not resolved.) The 
solid circles in Fig. 2 are the lineshape measured for an 
8CB + aerosil sample with ps = 0.025, the lowest aerosil 
density, at 287.8 K (T% A - 19.2 K). The larger width 
of the aerosil peak with respect to the resolution reflects 
the finite extent of the smectic correlations in the pres- 
ence of the aerosil gel, and the consequent destruction 
of the transition to a quasi-long-range ordered smectic 
state by the quenched disorder. The length scale for this 
short-range order depends strongly on ps, as revealed by 
the contrast between the low-temperature width for the 
ps = 0.025 sample and that for a ps — 0.220 sample, 
shown by the open circles in Fig. 2. 

To characterize these short-range correlations, we 
model the structure factor for the smectic ordering with 
the form 



S(q) 



g*mi) 

(1 + (<?|| " q m + qia + cqUV) 2 



(1) 



where the wave vector gii is along the smectic layer nor- 
mal and q± is perpendicular to it. The first term in S(q), 
an anisotropic Lorentzian with quartic corrections, has 
the same form as the structure factor in pure nematic 
8CB and represents the smectic susceptibility, with <ti 
equal to the magnitude of the susceptibility at the order- 
ing wavevector, qq. This term characterizes the critical 
thcraaal fluctuations on approaching the N-SmA transi- 
tioned. The second term, whose shape is given by the sus- 
ceptibility squared, is designed to account for the static 
fluctuations induced by the quenched disorder. This sec- 
ond term is motivated by studies on random field sys- 
tems. Such an expression, known as the disconnected 
susceptibility, has been shown to describe accurately the 
short-range correlations induced by_.the static fluctua- 
tions in random fLeid|Jsing magnctsEZI and has been jus- 
tified theoreticallyEa~E3 for these systems. A very similar 
expression has also been derived specifically fat. smectic 
liquid crystals in the presence of random fieldsES. In the 
limit of long-range order, the second term evolves into a 
Bragg peak, a 2 5(q\\ - q )5(q±]B. In applying Eq. (1) to 
the smectic correlations of 8CB in aerosil gels, £j_ and c 
are treated as functions of with and c(£n) set 

by their relations in pure 8CB, which are known to high 
precision^. 

The total measured scattering intensity is fit to the 
powder average of S(q) from Eq. (1), convolved with the 
resolution, plus a term to account for the background 
scattering from the aerosil gel shown in Fig. 1: 
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FIG. 3. Comparison of the smectic peaks measured for 
8CB in aerosil gels and those measured for pure 8CB above 
T% A = 306.97 K: (a) p s = 0.078 g/cm 3 (solid circles) and 
pure 8CB (open circles) each at 308.3 K; and (b) ps = 0.031 
g/cm 3 (solid circles) at 310.8 K and pure 8CB (open circles) 
at 310.5 K. 

/(g) = J dq' J dQS(ci')Res(q - q') + A(T)B(q) (2) 

The powder average, represented by the integral over 
solid angles, can be solved analytically for the S(q) form 
given by Eq. (1). The details of this calculation are pro- 
vided in the Appendix. The convolution with the reso- 
lution function, Res(q), was performed numerically. The 
shape of the background, B(q), was taken from measure- 
ments at a high temperature, where scattering from the 
aerosil gel dominates, and this quantity is multiplied by a 
temperature dependent parameter, A{T), which accounts 
for variations in the scattering contrast between the silica 
and 8CB due to differences in their thermal expansion. 
We found that A(T) typically differed from unity by less 
than 10% over the full temperature range of the measure- 
ments. 

At temperatures near and above T^ A , the measured 
lineshapes in the aerosil samples are described very well 
by the thermal fluctuation term alone (i.e., fitted a 2 are 
very small, and one can set a 2 = with almost no ef- 
fect on the other parameters). Thus, for T > T^ A , the 
smectic correlations in aerosil samples are very similar to 
the critical fluctuations of the pure liquid crystal. This 
strong similarity is demonstrated in Fig. 3 which com- 
pares the smectic peak for two values of ps with that 
measured for pure 8CB at temperatures several degrees 
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above T% A . 

We note that the lineshapes of the smectic peaks in 
Figs. 1-3 are nearly symmetric in q about the peak po- 
sition. Such highly symmetric shapes are, in general, 
not expected from powder averaging Eq. (1). We in- 
terpret this symmetry as a consequence of the relative 
values of £||, and c, particular to 8CB, which con- 
spire to produce such a shape. Indeed, this symmetry 
highlights the importance of including the quartic term 
(c > 0) in Eq. (1) to account correctly for the measured 
lineshapes. To illustrate this point, Fig. 4 shows the 
measured scattering peak from smectic fluctuations in 
pentylphcnylthiol-octyloxybenzoate (8S5) confined in an 
aerosil gel with ps = 0.030 at 2.56 K above the N-SmA 
transition temperature of pure 8S5. For a given the 
quartic term c is considerably smaller in 8S5 than in 8CB. 
This difference leads to a much more asymmetric shape 
to the powder average of the thermal term in Eq. (1). 
The solid line in Fig. 4 is a fit with Eqs. (1) and (2) with 
<i2 = Oo, so that only the thermal term is included in 
describing the smectic fluctuations above Tna- 

In measurements on pure 8CB above TjJ-j, we observe 
a very broad liquid structure peak near 0.3 A -1 . This 
feature vanishes in the nematic phase below T^j. Careful 
measurements also revealed evidence of this peak at high 
temperatures in 8CB confined in aerosil gels. We note 
that this feature introduces a small inconsistency in our 
treatment of the background in Eq. (2). Specifically, B{q) 
is obtained from the scattering intensity above Tjvj, but 
any high-temperature liquid crystal contribution should 
not be included as part of the background below Tni- 
However, we have found that the liquid crystal contribu- 
tions to B(q) are sufficiently small that including or ex- 
cluding them does not influence the characterization of 
the smectic peaks at q « 0.2 A -1 . The relative strengths 
of the aerosil and liquid crystal scattering can, however, 
provide a method for normalizing the scattering intensity 
and thus for comparing the intensities from different sam- 
ples. In particular, we decompose B(q) into two separate 
contributions: 

I(q, T > T° 7 ) = B(q) = v s B s (q) + v LC B LC (q) (3) 

where the shape of the liquid crystal contribution, 
Blc(q), is obtained from measurements on pure 8CB 
and the aerosil gel contribution, Bs{q), is assumed to 
have a simple power-law formEI. The ratio of the pref- 
&ctoTSJirc/vs scales with aerosil volume fraction as ex- 
pected^. The prefactor to the liquid crystal contribu- 
tion, vlc, which is proportional to the quantity of 8CB 
in the beam, provides a method for placing the scatter- 
ing intensities on an absolute scale. We define normalized 
strengths for the thermal and static fluctuation terms by 

erf = ax/VLC (4) 
a 2 = a 2 /v LC (5) 

which give the sizes of the two terms for each sample nor- 
malized by the scattering volume of 8CB for that sample. 
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FIG. 4. Smectic peak measured for 8S5 confined in an 
aerosil gel with p s = 0.030 g/cm 3 at 339.04 K, 2.56 K above 
the N-SmA transition temperature of pure 8S5. The solid line 
is a fit with Eqs. (1) and (2), using values of £_l(£m) and c(£n) 
from pure 8S5 and with a-z = 0, so that only the thermal fluc- 
tuations are included in describing the smectic correlations. 
The strongly asymmetric lineshape is a consequence of the 
powder averaging. 

Figure 5, which shows representative fits to the line- 
shapes of two samples with Eqs. (1) and (2), illustrates 
the good agreement we observe for this fitting form over 
the entire investigated range of densities and tempera- 
tures. The solid lines in the figure correspond to the full 
fit while the dotted, dash-dotted and dashed lines corre- 
spond to contributions made by each of the three terms - 
thermal fluctuations, static fluctuations, and background 
- respectively. As mentioned above, in the fits with 
Eq. (1), £j_ and c are treated as functions of £||, with 
an d c(£||) set by their relations in pure 8CB. Thus, 
each fit has only five free parameters - go, A(T), ctj., 
and 02 - with two parameters, £m and aifux, controlling 
the shape of the peaks. 

B. Smectic layer spacing 

The wavevector q in Eq. (1) characterizes the spac- 
ing of the smectic layers, and Fig. 6 displays the values 
of qo extracted from fits for three aerosil densities. At 
high temperatures, where the smectic peak has a large 
width and small amplitude, the scatter in qo is too large 
to determine confidently any systematic dependence on 
temperature or ps- At low temperatures, go appears to 
be virtually independent of ps and increases slightly with 
decreasing temperature. This temperature dependence is 
shown in Fig. 6 as a function of the difference in temper- 
ature from T* , the temperature below which the mag- 
nitude of static smectic fluctuations becomes non-zero. 
The procedure for obtaining T* and its value as a func- 
tion of ps are discussed below in Sec. HID. 

The temperature dependence for go was also measured 
for pure 8CB below T^ A . Due to small systematic un- 
certainties introduced by variations in the x-ray optics, 
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FIG. 5. Scattering intensity, I(q), for 8CB in confined 
aerosil gels, (a) p s = 0.025 g/cm 3 at 288.5 K, (b) p s = 0.025 
g/cm 3 at 307.6 K, (c) p s = 0.282 g/cm 3 at 300.6 K, and (d) 
ps = 0.282 g/cm 3 at 306.6 K. The solid line in each panel is 
a fit with Eqs. (1) and (2), and the dashed line is the back- 
ground contribution from the aerosil gel structure. In panels 
(a) and (c), the dash-dotted line is the contribution from the 
first term in Eq. (1), representing thermal fluctuations, plus 
the background, and the dotted line is the second term, repre- 
senting static fluctuations, plus the background. In panels (b) 
and (d), the temperatures are above a pseudotransition tem- 
perature T* defined in Sec. HID, and the smectic correlations 
are described entirely by thermal fluctuations. 
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FIG. 6. The ordering wavevector, go, extracted from the 
fits with Eqs. (1) and (2) for smectic correlations in 8CB con- 
fined in aerosil gels with three ps values as specified. Tem- 
perature is plotted as the difference from T* , the onset tem- 
perature for static fluctuations. Values for T* are given in the 
inset on Fig. 8. 

precise comparisons on an absolute scale between the val- 
ues of go for pure 8CB and those for the aerosil samples 
should not be made. However, q for the aerosil sam- 
ples at low temperatures and for pure 8CB in the smec- 
tic phase showed very similar temperature dependencies, 
and we estimate that for a given effective temperature 
difference, T — T*, the fractional difference in c/o between 
pure 8CB and any aerosil sample is not more than 0.2%. 
Thus, the presence of the aerosil has little, if any, effect 
on the partial bilayer smectic layer spacing in 8CB. 



C. Thermal fluctuations 

The amplitude of the thermal fluctuation term, erf, is 
displayed Fig. 7 as a function of temperature for several 
values of ps- Consistent with Fig. 3, the temperature 
dependence of erf at temperatures well above T* tracks 
that of the susceptibility for pure 8CB, which is shown 
by the solid line in the figure. However, unlike the pure 
system, where the susceptibility diverges at T^ A , erf for 
8CB confined in aerosil gels remains finite through the 
"transition" region. 

At temperatures below T* , erf exhibits a roughly 
temperature independent value. In Fig. 7, erf for the 
ps = 0.220 and ps — 0.282 samples are shown as free 
parameters to low temperatures in order to illustrate 
the very weak temperature dependence at low T. Be- 
low the pseudo-transition region near T*, the thermal 
term should be dominated by long wavelength excita- 
tions, and the scattering intensity from these contribu- 
tions should be dictated by the Bose occupation factor, 
which gives erf (T) oc T. Over the narrow absolute tem- 
perature range below T* (approximately 290 K to 307 
K), the roughly constant behavior for erf is consistent 
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FIG. 7. Amplitude of the thermal fluctuations, c>i , 
8CB confined in aerosil gels of several densities ps as spec- 
ified. Temperature is plotted as the difference from T* , the 
onset temperature for static fluctuations. Values for T* are 
given in the inset on Fig. 8. The solid line is the temper- 
ature dependence of the susceptibility for pure 8CB, with 
T*(pure) =T% A . 



with this expectation. On the basis of this observed low- 
temperature <Ji behavior in unconstrained fits, we have 
repeated the fitting with Eqs. (1) and (2) for each p$ for 
T < T* with <jf fixed at the average low temperature 
value for that density. This second iteration of fitting was 
designed to obtain better stability in the results for the 
other fit parameters (namely and £||) that are cou- 
pled to o-fS. The results for p s = 0.025,0.041,0.051, 
and 0.105 in Fig. 7 show erf 7 held constant at low tem- 
perature under this procedure. Near and above T*, we 



treat cr{ as a free, temperature-dependent parameter, 
leading to the behavior shown in Fig. 7. 
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D. Static fluctuations 
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FIG. 8. The integrated area, 0,% , of the static fluctuation 
term in S(q) for three values of ps'- 0.051 (solid circles), 0.105 
(open triangles), and 0.341 (solid squares). The values for 
ps = 0.051 g/cm 3 and ps = 0.105 g/cm 3 have been multiplied 
by 7 and 5, respectively, for clarity. The solid lines are fits 
with Eq. (6). The open circles in the inset are the values 
of T* , the onset temperature for static fluctuations extracted 
from these fits, as a function of ps- The solid circles in the 
inset are the peak positions in the specific heat as obtained 
from calorimetry studies [7]. 



quenched disorder. The open symbols in the inset of 
Fig. 8 are the resulting T* values as a function of ps- 
As the inset illustrates, the onset of static smectic order 
becomes increasingly suppressed in temperature as ps 
increases. This suppression in T* quantitatively tracks 
the decreasing position of the peak in the heat capacity 
with ps observed in calorimetry studiesQ and shown in 
the inset with solid symbols. Due to the considerable 
uncertainty in the values of T* extracted from the fits 
with Eq. (6), the nonmonotonic variation in "transition" 
temperatur&nwith ps that was observed calorimctrically 
for ps < 0.10 cannot be resolved. 



As mentioned above, the scattering intensity from 8CB 
in aerosil gels above T^ A is well described by thermal 
fluctuations alone, and the static term in Eq. (1) is effec- 
tively zero. With decreasing temperature, the contribu- 
tions from static fluctuations rise sharply from zero. Fig- 
ure 8 displays the temperature dependence of the static 
fluctuation term, a% , for three values of ps- The solid 
lines in the figure are the results of fits with the form 

a£ = B{T*-T)*, (6) 

where B, x, and T* depend on ps- Values for B, x, and 
T* are tabulated in Paper II. As the figure illustrates, 
this expression describes the temperature dependence of 
quite well. The significance of this functional form is 
discussed in Section IV. 

Fits with Eq. (6) provide unambiguously the onset 
temperature T* of static fluctuations induced by the 



E. Smectic correlation length 

Consistent with Fig. 3, the values of the smectic cor- 
relation lengths, £y , extracted from the fits with Eqs. (1) 
and (2) for T > T* are like those of pure 8CB and track 
its temperature dependence. Figure 9(a) displays £|| at 
temperatures above T*. Any possible systematic varia- 
tion in £|| with ps is overwhelmed by the scatter in the 
data. The solid line in the figure is £ii for pure 8CB@, 
where we use T^ A for T*. A detailed scaling analysis of 
the temperature dependence expected for ^ of 8CB + 
aerosils above T* is provided in Paper II. 

As with cf, £|| for aerosil samples fails to track the 
diverging behavior of the pure system through Tjva and 
remains finite to low temperature. At temperatures be- 
low T* , £|j approaches an essentially temperature inde- 
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FIG. 9. The smectic correlation length, £y , for 8CB con- 
fined in aerosil gels, (a) £y for T > T* compared with that of 
pure 8CB [solid line with T*(pure) = T^a] f° r values of ps as 
specified in Fig. 10. (b) £y values given to low temperatures 
for four values of ps- 



pendent value. Figure 9(b) shows the correlation lengths 
at low temperature for four densities of aerosil. The vari- 
ation of the low temperature correlation length, ^ T = 

£y(T << T*), with p s is described in Section IV. 

IV. DISCUSSION 

As the results outlined in Sec. Ill demonstrate, con- 
finement by aerosil gels strongly affects the smectic tran- 
sition in 8CB. The success of Eq. (1) indicates that these 
effects can be modeled as random fields. Consistent 
with theoretical predictions the random fields destroy the 
transition to a quasi-long-range-ordered smectic phase 
and replace it with the formation of short-ranged cor- 
relations. In this Section, we investigate the relationship 
between the strength of this disorder and the correspond- 
ing short-ranged smectic order, including possible con- 
nections with theoretical predictions, as well as aspects 
of the pseudocritical behavior exhibited by 8CB confined 
in aerosil gels. 
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FIG. 10. Amplitude of the thermal fluctuations, erf, ver- 
sus correlation length £y for 8CB confined in aerosil gels. The 
plotted values span the full experimental range of ps and tem- 
perature. The solid line is the relation erf oc 



A. Pseudocritical behavior 

As Figs. 7 and 9(a) illustrate, the smectic scattering 
in 8CB confined in aerosil gels above T* resembles, to 
a first approximation, thermal critical fluctuations unaf- 
fected by quenched disorder. This close connection be- 
tween the thermal term in Eq. (1) for 8CB + aerosils and 
the corresponding term for pure liquid crystals in the ne- 
matic phase is strengthened by considering the scaling 
relation between erf and In a pure isotropic critical 
system for T > T c , one has o\ oc t~ 7 and £ cx t~ v which 
leads to o\ oc = (, 2 ~ v . For liquid crystals, which 
exhibit critical anisotropy in the behavior of the correla- 
tion lengths, one obtains <j\ oc £?/ V ^ where 7/V|| = 2 — jjy. 
This scaling relation between the thermal susceptibility 
and the parallel correlation length for 8CB + aerosils can 
be written as 



(7) 



where the subscript on r\ has been dropped for conve- 
nience. Figure 10 displays a log-log plot of erf versus 
£jj for a series of samples of varying ps and a line rep- 
resenting <7i oc £|. As the fi gure demonstrates, the scal- 
ing relation holds for 8CB confined in aerosil gels over 
the entire experimental temperature range, both in the 
pseudo-critical regime above T* and at low temperature 
where both erf and £|| have only weak temperature de- 
pendence. Further, the scaling amplitude, A, shows no 
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systematic variation with ps- Indeed, the successful col- 
lapse of data from different samples in Fig. 10 suggests 
the validity of the normalization procedure for using 
vlc- Also clear from the figure is the fact that 2 — rj is 
close to 2 for 8CB + acrosils. The mean field and the 
Gaussian tricritical values of 2 — rj are 2.00, and the 3D- 
XY value is LQ62. Furthermore, the value of 7/^1 is 1.88 
for pure 8CBc3. The correspondence between the scaling 
given in Eq. (7) applied to 8CB + aerosil gels of varying 
ps and that for pure liquid crystals is discussed below. 

This pseudocritical behavior of 8CB confined in aerosil 
gels also extends below T*, where c in Eq. (1) approaches 
a small constant value, and 0,2 becomes proportional to 
the integrated intensity of the static fluctuation term 
in I(q). The growth of this intensity, shown in Fig. 8, 
strongly resembles that of an order parameter squared, 
as the success of the fits of a 2 with Eq. (6) illustrates. 
Figure 11(a) shows the values of the effective exponent 
x extracted from these fits as a function of ps- For pure 
liquid crystals, the critical behavior observed at the N- 
SmA transition is affected by the proximity of the higher 
temperature nematic to isotropic transition through a 
coupling between the nematic and smectic order-param- 
eters (e.g. de Gennes coupling; see Paper II)E3. The 
strength of this coupling depends on the magnitude of 
the nematic susceptibility at Tna, which for pure liquid 
crystals depends on the width of the nematic range and 
can be roughly parameterized by the MacMillan ratio, 
Rm = Tna/Tni. When Rm becomes small (Rm ~ 0.7), 
the N-SmA transition is observed to approach 3D-XY 
behavior but still exhibits a small anisotropy in the criti- 
cal correlationsEJ, while it approaches tricritical behavior 
and then becomes first order as Rm — > 1» Values of 
2/3, determined through the Rushbrooke scaling equality, 
2/5 = 2 — a — 7EZI, for a series of pure liquid crystalsEa 
with varying Rm are included in Fig. 11(a). The effective 
8CB + aerosil exponents, x, span the range from near the 
tricritical value 2/3 = 0.5 (like pure 8CB) at small ps to 
3D XY-like, 2(3 = 0.691, at large p s . With increasing 
ps, the quenched disorder has the apparent effect of sup- 
pressing the nematic susceptibility and thus the coupling 
that drives the critical behavior of the N-SmA transition 
in pure liquid crystals away from 3D-XY universality. 

This trend agrees very closely with the corresponding 
behavior observed in specific heat studies of 8CB with 
dispersed aerosilal Figure 11(b) shows the values of the 
effective heat capacity exponent a from the calorimetric 
study versus ps- As the figure illustrates, a decreases to- 
ward the 3D-XY value with increasing ps- The solid cir- 
cles in Fig. 11(b) are the critical exponents measured for 
pure liquid crystals as a function of Rm- The same lin- 
ear scaling between ps and Rm successfully collapses the 
aerosil data onto the pure liquid crystals in Figs. 11(a) to 
11(b) and suggests the notion of an effective MacMillan 
ratio for 8CB in aerosil gels, 
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FIG. 11. (a) The effective exponent x from fits to 012 (T) 
with Eq. (6) for 8CB+aerosil samples (open circles) versus 
ps. Also given are values of the squared order parameter ex- 
ponent, 2/3 (filled circles), for a series of pure liquid crystals 
with varying MacMillan ratios, Rm = Tna/Tni, as deter- 
mined from data given in Ref. [34]. The effective exponent 
x spans the range from near a tricritical value, 2/3 = 0.5, at 
small ps to a 3D XY-like value, 2/3 = 0.691, at large ps. (b) 
the effective exponent, a, from calorimetric studies on 8CB 
confined in aerosil gels [7] (open circles) along with values of 
a (filled circles) for a series of pure liquid crystals with vary- 
ing Rm, as given in Ref. [34]. (c) the effective exponent 2 — 77 
obtained from fits with the form — A(?, ~ v versus ps (open 
circles) and the critical exponent j/vu (=2 — rj«) for a series 
of pure liquid crystals as a function of Rm- Note that the 
linear scaling between ps and Rm, given by Eq. (8), is the 
same for all three plots. 
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This general relationship between R e // (sil) and ps for 
two distinct pseudo-critical behaviors strengthens the 
idea that the presence of quenched disorder suppresses 
the nematic susceptibility in 8CB near T*. This sup- 
pression likely results from the role of surface anchor- 
ing on the gel strands in changing the properties of the 
nematic since the observed_ratio Tna/Tni for 8CB + 
aerosils varies little with psu- Figure 11(c) displays val- 
ues of 2 — m extracted from fits with Eq. (7) to erf and 
£|| values obtained for different ps samples. Also shown 
in Fig. 11(c) are values of (= 2 — 7711) for pure liquid 
crystals as a function of Rm- Due to the small difference 
in 2 — 7711 between the tricritical and 3D-XY values and 
the large scatter in the results both for pure liquid crys- 
tals and for 8CB + aerosil samples, systematic trends 
like those seen in Fig. 11(a) and 11(b) are less clear for 
2 — rj. Tabulated values for 2 — rj are given in Paper II. 



B. Low temperature correlations 

As Fig. 9(b) demonstrates, the correlation length for 
smectic order in 8CB confined in aerosil gels saturates at 
low-temperature. Figure 12(a) shows these low tempera- 
ture parallel correlation lengths, ^ T , as a function of p$- 
The solid line in Fig. 12(a) shows the trend in the nematic 
correlation lengths of the homolog 6CB with dispersed 
aerosils, as measured by Bellini et alS3 with static light 
scattering. These light scattering studies demonstrate 
that the nematic order in a liquid crystal with dispersed 
aerosil breaks up into very large, but finite, domains. 
Since smectic layering forms within regions with nematic 
order, the nematic domain sizes set upper limits for the 
range of the smectic correlations. However, remains 
well below these limits for all ps , demonstrating the sen- 
sitivity of the smectic phase to quenched disorder. The 
dashed line in Fig. 12CaX.is the mean aerosil void size, 
Iq, as a function of psQ'Ex Over the range of densities 
studied, £,^ T exceeds Iq by a factor of three or more, con- 
sistent with an intermediate strength disordering field. 
Paper II provides a detailed discussion of the role of 1$ 
in setting the range of the smectic correlations. 

Figure 12(b) displays af /cf as a function of ps for 
T = T* - 6 K. We plot this ratio rather than erf or af 
alone to characterize the magnitude of the fluctuations in 
order to remove any errors that the uncertainties in vlc 
might introduce through the normalization. Because the 
exponent, x, characterizing the growth of af below T* is 
only weakly dependent on ps, the trends observed for this 
ratio are insensitive (on a logarithmic scale) to the tem- 



perature difference T* — T used to evaluate 



Figure 



12(c) shows the quantity /erf, where £ = (£||£i) 1/3 
is the mean correlation length at low temperature. This 
quantity, which is approximately the ratio of the inte- 
grated intensities of the thermal and disorder terms at 
low temperature, shows little variation with ps- 
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FIG. 12. (a) Low temperature correlation length, ^ T , as a 
function of ps- The solid line is the nematic correlation length 
(for 6CB) measured by static light scattering [10], and the 
dashed line is the mean aerosil void size, lo- The long-dashed 
line represents the form oc p^ 1 . (b) The ratio af /erf (low 
T) at the fixed temperature difference T* — T = 6 K as a func- 
tion of ps- The long-dashed line is the relation af /af (low 
T)oc p%. (c) The ratio af £ 3 /erf (low T) at the fixed temper- 
ature difference T* — T — 6 K as a function of ps- 
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C. Comparisons with theory 

As mentioned in the introduction, random field mod- 
els have served as important examples for understanding 
possible consequences of quenched disorder, and exten- 
sive theoretical work has been devoted to random field 
systems. For the case of transitions that break a contin- 
uous symmetry in three dimensions, like the nematic to 
smectic-A transition, the domain wall energy arguments 
of Imry and Ma conclude that the-transition is unstable 
to arbitrarily weak random fieldsEj. Subsequent theo- 
retical work has rigorously established this conclusion^, 
consistent with the short-ranged smectic order measured 
for all 8CB + aerosil samples well below T^ A . For quan- 
titative comparisons with these short-range smectic cor- 
relations two theoretical efforts provide guidance. In the 
first, Aharony and Pyttai-3 apply scaling arguments to 
predict the low temperature form of the structure factor 
for systems with transitions breaking continuous sym- 
metries in the presence ©f random fields. In the second, 
Radzihovsky and TonertH develop a detailed theory for 
smectics with random quenched disorder. 

The strength of the disorder imposed on smectic 8CB 
by confinement in aerosil gels clearly correlates with the 
density of the aerosil; however, the precise relation be- 
tween density and disorder is unknown. Within the sim- 
plest assumption, confinement in the gels induces random 
fields whose disorder strength is proportional to the gel 
density. With this assumption, the scaling arguments of 
Aharony and Pytte predict £ oc A -1 cx p^ 1 , where A is 
the strength of the random field disorder. The exponent 
— 1 is obtained from —l/(d c — d) where d = 3 is the spa- 
tial dimension and d c = 4 is thc-iowcr critical dimension 
for the random field XY modeH. The long-dashed line 
in Fig. 12(a) is the power-law form oc pg 1 . Over 
the range of ps covered experimentally, measured values 
of appear roughly to follow this behavior. A best 
fit to a power-law form over the full range of ps gives 
cx Pg 1 ' 1 ^ ' 2 - The scaling arguments of Aharony and 

Pytte further lead to /cr^ oc A 3 . The long-dashed line 
in Fig. 12(b) is this relation, again assuming A oc ps- As 
with £||, ci2 /<J\ follows a trend crudely consistent with 
the predicted power-law behavior, but clear deviations 
from the predicted form are also apparent. One possi- 
ble source for these discrepancies could be a more subtle 
relation between ps and disorder strength. Paper II elab- 
orates on this issue. 

In order to make comparisons that are free of the de- 
pendence of disorder on ps, we plot in Fig. 13 /oi 
against £. In this plot we use £ rather than £y since 

3 

the scaling theory is for an isotropic system, and £ rep- 
resents the correlation volume. According to the scal- 
ing theory, /cr™ should vary with correlation length 

3 

like £ , while a best fit shown by the solid line gives 

-r-2.68±0.26 , • , , , 

£ . 1 his relatively good agreement with the the- 

ory is also reflected in the weak variation with disorder 




° log 10 (S) 

FIG. 13. The ratio a^/crf' (low T) at the fixed temperature 
difference T* — T — 6 K as a function of £. The solid line is 
a power law fit giving exponent value of -2.68. 

strength of a^tt/ai shown in Fig. 12(c). 

A second possible source of discrepancy between the 
smectic correlations in 8CB + aerosil and the scaling ar- 
guments of Aharony and Pytte is that additional disor- 
dering effects enter the confined smectic beyond those 
of the random field 3D XY system. In their discussion 
of smectic liquid crystals with quenched disorder, Radz- 
ihovsky and Toner identify two sources of disorder, tilt 
disorder that couples to orientation (ie, the nematic di- 
rector) and a layer displacement disorder that couples 
directly to the smectic order parameter. While the lat- 
ter corresponds to the disorder of the XY model, the 
former is predicted to be the dominant disorder for the 
smectic. An important prediction of the work of Radzi- 
hovsky and Toner is that the disordered smectic should 
possess anomalous elasticity, and one consequence of this 
prediction is that the correlation length ratio, is 
reduced from that of the pure system. In our modeling of 
the smectic correlations with Eq. (1), we find that main- 
taining equal to its values in pure 8CB produces 
good fits to the data. However, the powder nature of the 
measured lineshapes creates considerable uncertainties in 
so that we can not exclude some reduction in the 
correlation length ratio. Despite this inconsistency, the 
prediction from Radzihovsky and Toner that oc A -1 
is roughly consistent with our results in Fig. 12(a), as- 
suming again A oc ps holds approximately. 

Radzihovsky and Toner further predict that £(T) oc 
BiT) x / T at low temperature, where B{T) is the layer 
compression modulus of the pure liquid crystal. Within 
this theory, 1/r > 1/2 is required for a smectic Bragg 
glass phase to be stable. Published data for B{T) of 
8CB extend only 3 K below Tjy^Ea, so a direct test of 
this scaling prediction is difficult. However, the essen- 
tially temperature independent behavior of £n at low tem- 
perature shown in Fig. 9(b) indicates that such scaling 
does not hold for 8CB in aerosils. We note that this 
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conclusion differs from that reported for 8CB confined 
in aerogels, in which the analysis of those x-ray scatter- 
ing results yielded £(T) consistent with£(T) oc B^) 1 / 7 
and with marginal Bragg glass stability!! This contrast 
between the aerosil and aerogel results is somewhat sur- 
prising since the validity of the prediction of a Bragg 
glass is within a regime of weak disorder, and the aerosil 
gels impose a more gentle perturbation on the smectic 
than aerogels. One source for the different £\\{T) behav- 
iors could be the partial compliance of the aerosil gel to 
the elasticity of the liquid crystal, a feature which is not 
incorporated into the theory and which may complicate 
comparisons with it. However, a number of differences 
also exist between the approach used .for analyzing the 
smectic lineshapes in the aerogel studycrtil and the analy- 
sis presented here for aerosil samples. For 8CB + aerosil 
data, analysis like that used in the aerogel study has 
strong consequences for £ (T) . Further study of other liq- 
uid crystals would clearly help in determining whether 
the contrast between the aerogel and aerosil systems in 
their agreement with Bragg glass stability is a conse- 
quence of elastic coupling that is not considered in the 
theory or, perhaps, an artifact of the differing approaches 
to the lineshape analysis. Such studies would also be in- 
valuable in testing the generality of trends observed for 
8CB, such as the crossover to 3D XY behavior shown 
in Fig. 11, with the introduction of quenched disorder 
through confinement in aerosil. 
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APPENDIX: ANALYTIC POWDER AVERAGING 
OF THE SMECTIC STRUCTURE FACTOR 

The x-ray structure factor for the smectic fluctuations 
in 8CB confined by aerosil gels, given in Eq. (1), in- 
cludes two terms. The first term (LFC) is an anisotropic 
Lorentzian with a 4th order correction in the transverse 
direction that accounts for critical thermal fluctuations. 
The second term (LFC 2 ) is proportional to the ther- 
mal term squared and accounts for random-field induced 
static fluctuations. Although an anisotropic Lorentzian 
raised to a power along the transverse direction, S(q) = 

a/{q(q ll -q ) 2 + (l + e±ql) 1 - ,1±/2 }, is equally sue- 
cessful in describing the thermal fluctuations^, we use 
the LFC because its powder average can be calculated 



analytically. The structure factor for a powdered sam- 
ple is equivalent to that for a single-domain sample aver- 
aged with equal weight over all orientations of the smectic 
wave vector, qo, relative to the scattering wave vector, 

q, 

SP o wd er {q) = ^_J 5(q)< (Q) 

For the LFC term, the integral has the form 



T 



(10) 



(A + 2Bf, + C> 2 + DM 4 ) ' 



where 



A = 1 + rifgg + + c&q 4 , 



-II 

B = -£ 2 q q, 

*2 f2 \ 2 



(11) 



c = iq-&)(?-2c& q *, 

D = c£l<? 4 . 
This integral in turn can be written as 

sss-h) -5 />n <^)p <«> 

where ot^ are the roots with multiplicity m of the quartic 
equation: 



A + 2B(i + C/i 2 + = 0. 



(13) 



With the integrand in this form, it can be expanded into 
a sum of terms, 
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where the coefficients Cy are obtained through the ex- 
pansion process and are listed in Table I. At this stage, 
each term can be evaluated analytically using elementary 
definite integrals, such as 
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(15) 



where a = an + iai. 

The analytic powder averaging of the LFC 2 term fol- 
lows the same approach: 

s powd "'{q) — -i- f dn 

LFC 2 4-ir J 



{i+ 5 2 ( 9|| - 50 ) 2 + £M + <4}' 

*/>n,^ as) 
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As with the LFC term, this integral can be expanded into 
a series of elementary definite integrals. The coefficients 
for the expansion of the LFC 2 integral are listed in Table 
II. 
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TABLE I. Coefficients for terms potentially appearing in the expansion of the integral for the powder average of the LFC. 
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TABLE II. Coefficients for terms potentially appearing in the expansion of the integral for the powder average of the LFC . 
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